GLOBAL DYNAMICS BELOW THE GROUND STATE ENERGY 
FOR THE ZAKHAROV SYSTEM IN THE 3D RADIAL CASE 
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Abstract. Wc consider the global dynamics below the ground state energy for 
the Zakharov system in the 3D radial case. Wc obtain dichotomy between the 
scattering and the growup. 
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1. Introduction 

In this paper, we continue our study [7] on the global Cauchy problem for the 3D 
Zakharov system 

iu — Au = nu, 
h/a 2 — An = — A\u\ 2 

with the initial data 

u(0, x) = uo, n(0, x) = no, n(0,x) = ni, (1.2) 
where (u, n)(t,x) : R 1+3 — > C x E, and a > denotes the ion sound speed. It 
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preserves ||it(i)||x,a and the energy 



E = I \Vu\ 2 + \^HIM±K _ n \ u?dx , (1 . 3) 



where D := v— A, as well as the radial symmetry. 
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This system ( 11. ip in d dimensions was introduced by Zakharov [2T| as a mathe- 
matical model for the Langmuir turbulence in unmagnetized ionized plasma. It has 
been extensively studied. Local wellposedness (without symmetry) is well known. 
For example, the well-posedness in the energy space was proved in [I] for d = 2, 3 
and in [6] for d = 1 , and in weighted Sobolev space in [11] . It has been improved to 
the critical regularity in [6j |3] for d = 1,2, and to the full sub critical regularity in 
[6l[2]for<i>4,<i = 3. The well-posedness for the system on the torus was studied 
in [201 [T2"] . These results except for [TT] follow from the iteration argument using 
Bourgain space, where the estimates depend on a, while in [TT] the well-posedness 
is obtained uniformly for a. For more results on the subsonic limit to NLS (as 
a — > oo), see [J8J [TBI H3]- Concerning the long-time behavior, Merle [JT] obtained 
blow-up in finite or infinite time for negative energy (which we will call grow-up for 
brevity), while the scattering theory was studied in [T9], [5], [TT] , dealing with solutions 
for given asymptotic free profiles. Recently, in [7] the authors obtained scattering 
for radial initial data with small energy in the 3D case, by using the normal form 
reduction and radial-improved Strichartz estimates. The purpose of this paper is to 
consider the global dynamics for larger data under the radial symmetry. 

To simplify the presentation, we rewrite the system into the first order as usual. 
Let N := n — iD~ 1 h/a. Then (11.11) can be rewritten as 

(id t - A)u = (UN)u, [id t + aD)N = aD\u\ 2 , (1.4) 

with initial data (uq, N ) G H 1 x L 2 . It has the conserved mass 



M(u) : = 




and the Hamiltonian 

Ez{u, N) := I Vf + \f-^ dx = Es{ u ) + llN-luny4, (1.6) 
7r3 I 4 I 

where Es(u) denotes the Hamiltonian for the cubic NLS (the limit a — » oo) 

(id t - A)u = \u\ 2 u, (1.7) 

namely 

E s (u) := [ ^ - ^dx. (1.8) 

Let Q be the ground state for NLS (11.71) . that is the unique positive radial solution 
for the following equation 

-AQ + Q = Q 3 , (1.9) 

which minimizes the action 

J(Q) :=E S (Q) + M(Q) (1.10) 

among all nontrivial solutions of (jl.9p (see, e.g., [8] for further properties of Q). For 
A > 0, let 



Q\(x) := XQ(Xx), 



(1.11) 
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then we have 

-AQ X + \ 2 Q X = Ql M(Q X ) = X~ 1 M(Q), E S (Q X ) = XE S (Q). (1.12) 

Thus the Zakharov system (11.41) has the following family of radial standing waves 

(u,N) = (e^- x ^Q x ,Ql), (1.13) 

where A > and AgR can be chosen arbitrarily. 

The goal of this study is to determine global dynamics of all the radial solu- 
tions "below" the above family of special solutions, in the spirit of Kenig-Merle 
|10j . namely the variational dichotomy into the scattering solutions and the blowup 
solutions. Such a result has been obtained for the limit equation (jl.7p by Holmer- 
Roudenko j8] in the radial well as in the nonradial case [9]. For the di- 

chotomy, we need to introduce another functional (for NLS), which is the scaling 
derivative of the action J: 

K(ip) := d x \ x=l J{\ d / 2 ^{\x)) = [ \V V \ 2 -^dx. (1.14) 

We would like to get the same result as in [8] for NLS, but by the virial argument as 
in [TJ] we can only prove grow-up, due to the poor control of the wave component 
N. In fact, existence of any blowup in finite time is still an open question for the 
3D Zakharov system. The main result of this paper is 

Theorem 1.1. Assume that (w , N ) G if 1 (lR 3 ) x L 2 (M 3 ) is radial and satisfies 

E z (u , N )M(u Q ) < E S (Q)M(Q). (1.15) 

Then we have 

(a) if K(v,q) > 0, then (jl.4p has a unique global solution (u,N), which scatters 
both as t —7- oo and as t : — ^ — oo in the energy space. More precisely, there are 
(u±,N±) G H 1 x L 2 such that 

\\(u(t),N(t))-(e- UA u ± ,e UaD N ± )\\ mxL2 ->0 (t ±oo). (1.16) 

(b) if K(u$) < 0, then (11.41) blows up in either finite or infinite time, in the sense 
that sup 0<t<r » ||(«, A0||Hi(R3) xi 2 (]R 3) = oo = sup r>t<t<0 ||(n,iV)|| H i (IR 3 )xL 2 (R 3 ) , where 
(T*,T*) is the maximal interval of existence. 

Remark 1.1. 1) Assuming K(u ) = and (11.151) . one can actually get by variational 
estimates that u = 0, so u = and iV = e ltaD N , see Section 2. 

2) The condition (I1.15P is sharp in view of the standing wave solutions (I1.13p . 

The difficulty for the scattering even for small data can be observed by comparing 
the time decay with the NLS of general power nonlinearity 

iu-Au=\u\ p u, u:R 1+d ^C. (1.17) 

It is well known that the scattering for NLS requires p > 2/d, corresponding to the 
time integrability of the optimal decay of to the potential 

|||nn| Li o~|t|-*/ 2 , (1.18) 
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while the scattering in H s for any s has been proven only for p > 4/d. The 3D 
Zakharov system would be on the boarderline in the above sense, since the potential 
n can decay only by 

IMU? ~ (i.i9) 

as it is solving the 3D wave equation. This suggests that the decay estimates are 
far insufficient for the scattering in H 1 , and so it is essential to exploit nonlinear 
oscillations, e.g. by the normal form. This part for small radial data has been 
resolved in the previous paper [7]. Hence our main task in this paper is to carry 
out the Kenig-Merle approach [10] in accordance with the normal form. Since the 
normal form produces nonlinear terms without time integration, we need to modify 
Kenig-Merle's formulation, as well as some estimates in [7J. As a crucial ingredient 
for that approach, we will derive a virial identity, which is slightly different from 
Merle's one in [H] and more suitable for the scattering. 



2. Hamiltonian and variational structures 

2.1. Virial identity. We derive a virial identity on M. d , which is slightly different 
from [13]. Recall that the Zakharov system can be rewritten in the Hamiltonian 
form 

dt(fy =JE' z (u,N), (2.1) 

where J and E' z denote the symplectic operator and the Frechet derivative given by 

i \ _ (E' s (u) - (m - \u\ 2 )u\ _ ( -{m)u 

2iaDJ> ^z{u,!V)-y (iV-|n| 2 )/2 J ~ \(N - \u\ 2 )/2 

Let A be the generator for the family of scaling transform^ 



hence we have 

.4 



A* = - 



(2.3) 



x-W + d/2 

x- V + (d+ l)/2 y 

x-V + d/2 

x-V + {d-l)/2 / 

Let w := (u,N), J := J -1 and denote the real part of L 2 inner product by (-|-). 
Then the virial identity for the Zakharov system is given by 

d t (Jv\Av) = {3v\Av} + (Jv\Av) = (v\(3*A + A*J)v) 

= 2{v\ J*Av) = 2(JE' z (v)\3*Av) = 2{E' z {v)\Av) 

= 2d x=1 E z (S x v) = 2d x=1 [E s (S x u) + WX^N - X d/2 \u\ 2 \\l/^] (2 " 4) 
= 2K («) + \\\N - \u\ 2 \\l - d -^{N - H 2 |H 2 ), 



1 The order of scaling, i.e. the exponents | and ^Jl, is the unique choice such that (|2.5j) holds. 
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where we used for the third equality that 

-f* A * _ ■ fx ■ V + d/2 | — /It (n k\ 

JA - Z y o 2a(x-V + (d+l)/2)Dj ~ AJ - ^ 

Therefore, we have proved 

Lemma 2.1 (Virial identity). Assume v = (u,N) is a smooth decaying solution to 
Zakharov system (11. ip . Then 

d t (3v\Av) = d t \(u\ird r u) + — (N\ird r D~ l N)] 

= 2K{u) + i||Af - |«| 2 ||I - — y— - \u?\H)- 



The virial identity by Merle [M] is slightly different from the above one. In our 
notation, it can be written as 

d t \(u\ird r u) - -(3lN\rd r D- 1 $SN)] 

= 2K{u) + ^\\N- \u\\\l - (d-l)\m\\l ( 2 - 7 ) 
= 2dE z (u,N) - (d-2)\\Vu\\l- (d- l)piV||2. 

The left hand side differs from ( 12. 6 p since ird r D~ l is not self-adjoint, but ir(d r + 
(d — 1)/2)D~ 1 is so. Precisely, the difference is 

d t ^l(m\D-^N) = i^(iV - M 2 |iV> + (d- (2.8) 

let 2 
The advantage of our identity is that it is monotone both in the scattering region 
[K > 0) and in the blow-up region (K < 0), as we will show in the next section, 
while ( 12. 7p is not monotone when u(t) and n(t) are very small compared with n(t). 
Although Merle's identity is more convenient in the blow-up region, our identity can 
also be used there, as we will see in Section |3j 

2.2. Variational estimates. In the 3D case d = 3, the cubic nonlinearity is L 2 - 
supercritical and H 1 subcritical. Hence Q is obtained by the constrained minimiza- 
tion 

J(Q) = inf{ J{<p) | ^ <p, K{ip) = 0}. (2.9) 

Indeed, Q is the unique minimizer modulo the phase e ld and spatial tranlation. By 
scaling, we also have for any A > 

XJ(Q) = J X (Q X ) = inf{ J A (y?) | ^ tp, K{tp) = 0}, J x := E s + A 2 M, (2.10) 

and Q\ is the unique minimizer modulo phase and translation. 

Lemma 2.2. Assume that (u,N) is a solution to (II. 4p with maximal interval I 
satisfying 

E z (u,N)M(u) < E S (Q)M(Q). (2.11) 

Then for some A > we have E z (u, N) + X 2 M{u) < XJ(Q). Moreover, either u = 
on I, or K{u{t)) ^ for all tel. In other words, K{u{t)) does not change its sign 
on I . 
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Proof. From (12.91) . we have J(Q) = infA>o J(Q\), and thus d\\\=iJ(Q\) = 0. This 
implies 

J 2 (Q)/A = E S (Q)M(Q). 
Thus we see that there exists A > such that 

E z (u, N) + \ 2 M{u) < J X {Q X ) = XJ(Q). (2.12) 

Since J\(u) < Ez(u,N) + X 2 M(u), by the variational characterization of Q\, we 
have at each t G /, 

K(u(t)) = ^ u(t) = 0. (2.13) 
If K(u(t )) = for some t G /, by uniqueness we have u = 0. □ 

Corollary 2.3. Assume that (u, N) is a solution to (11.41) with maximal interval I 
satisfying for some A > 

E z (u,N) + \ 2 M{u) < XJ(Q), K(u ) > 0. (2.14) 

Then I = (— oo, oo) ; and moreover, 

E z (u, N) + \ 2 M(u) ~ \\u\\ 2 m + \\N\\ 2 L2 ~ \\uo\\ 2 H i + ||iV ||£ 2 . (2.15) 

where the implicit constant depends only on A and J(Q) . 

Proof. From Lemma 12.21 (b) we get that if K(uq) = 0, then u = 0, and hence this 
case is trivial. Thus we may assume K(uq) > 0, hence K(u(t)) > by Lemma [2.21 
(b). From the assumption, we get (I2.15P immediately from 

XJ(Q) >E z (u, N) + \ 2 M{u) - K(u(t))/3 
i \ 2 i 

IIV7 1 1 2 i A II ||2 , 1 II at | 1 2 1| 2 
= gl|Vw|| 2 + y||u|||+-||A^- \U\ \\ 2 , 

and the Sobolev inequality ||n|| L 4<||n||^i. So (u,N)(t) is a priori bounded in H 1 xL 2 , 
and thus by the local wellposedness we have I = (— oo, oo). □ 

So far, the global well-posedness of part (a) of Theorem 11.11 is proved. It remains 
to prove the scattering and part (b). For both purposes, the virial estimates play 
crucial roles. Unlike the NLS case, it is not at all obvious that virial for (II. 4p is 
monotone. The following lemma is our key observation 

Lemma 2.4. Let y? G H 1 ^ 3 ), A > and v > satisfy 

E s {y) + X 2 M( V ) + ^ < J A (Q A ). (2.16) 



Then we have 

(K(<p) > =^ AK(<f) + v 2 > VEu\\<p\\l 
\K(ip) < =^ AK(cp) + v 2 < -2v\\<p\\l 



(2.17) 



Proof. First, if K(tp) = then v = and the conclusion is trivial. Hence we may 
assume K(ip) ^ as well as v > 0. Next by the scaling (if, v) (Xif(Xx), \/~Xv), we 
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may remove A or assume A = 1. Then the energy constraint becomes J(ip) + u 2 /A < 
J(Q). Now consider the L 2 scaling of (p, S^Lp = fx d / 2 ip(fix) and 



A^) = ^liv^ + J|MI1-^IMI1 



2 " rnz 2" r " z 

3 

~4 



(2.18) 



There is a unique < \i ^ 1 such that 

l|V^ = ^|M|t (2.19) 

which is equivalent to K^S^if) = 0. Then the variational characterization of Q 
implies J(S fl <p) > J(Q), and so 

£ < J(^) - JM = ^liv^ - £^±\\<p\\i, 

(^-1)^ + 2) (2 - 20) 

= § IML, 

where (12 . 19[) is used in the last step. Let X := \\(p\\ 2 /u. Then the above inequality 
is rewritten as 

\fi- l\y/fi + 2X > V2. (2.21) 
Hence it suffices to estimate, under the above constraint, 

AK{ ^ t ^ = 3(/i - l)X + 1/X =: /(X, p). (2.22) 

For K(f) > 0, or equivalently fi > 1, f(X,fj) is increasing in X unless 



<-^r\l^> (2-23) 



3(/i — 1) ii — 1 V A* + 2 
which is solved > (a/33 — l)/2. In the latter case, we have 



3(/x - 1)X + 1/X > 2v/3(/x - 1)X/X > \/6, (2.24) 
since /i > 3/2. Otherwise, the minimum is attained at the boundary and equal to 



+ " = : (2 - 25) 



which is increasing^ in /i > 0, hence 6(/i) > 6(1) = a/6. 

For K(tp) < 0, or equivalently < ji < 1, —f(X,/i) is increasing in X, so its 
minimum is attained at the boundary and equals to 



-/(- J^ro./i) = %) > 6(0) = 2. (2.26) 

1 - /i V H + 2 

Therefore, the proof of the lemma is completed. □ 



2 This can be checked by computing d ^ ^ . 
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Remark 2.1. Applying the lemma above by letting 

v:= \\N- \u\ 2 \\ 2 , (2.27) 

we get from Lemma I2TT1 that the virial (3v\Av) is monotone in our consideration. 
This fact will play crucial role in our consequent analysis. 

3. Growup at infinity 

This section is devoted to prove part (b) of Theorem 11.11 We assume that under 
the assumption of part (b), the solution exists for all t > 0. We will show that 

SUP \\(U,N)\\ H 1 {R 3 )><L 2 {R 3) = oo. 
t>0 

3.1. Localized virial. Let X = X* be the operator of smooth trancation to \x\ < R 
by multiplication with ipn(x) = ifj(x/R), where -0 G C^°(]R 3 ) is a fixed radial function 
satisfying < ij) < 1, d r ip < 0, ip(x) = 1 for \x\ < 1 and ip(x) = for \x\ > 2. We 
consider the localized virial quantity in the form 

V R (t) := (Jv\(AX + XA)v). (3.1) 
Then similarly to the non-localized virial identity, we can compute 

V R = (E' z (v)\(AX + XA + AJXJ + J*XJ*A)v). (3.2) 
Putting v := n — \u\ 2 , the right hand side can be written componentwise 
V R =(E' s (u) - uu\2A Xu + 2XA u) 

+ (u/2\(XA 1 + A ± X + DXD^A^ A 1 DXD~ 1 )(u+ \u\ 2 )), 
where Aj := x ■ V + (d + j)/2. The right hand side is decomposed into the NLS part: 

NS := (E' s {u)\2A Xu + 2XA u), (3.4) 
the quadratic terms in v: 

QN := (u/2\(XAx + A X X + DXD~ X A X + A x DXD~ x )v), (3.5) 

and the cubic cross terms: 

CC :=(-uu\2A Xu + 2XA u) 

+ (u/2\(XA 1 + A X X + DXD~ 1 A 1 + A^XD' 1 )^ 2 ) , 



(3.3) 



(3.6) 



i.e., Vr = NS + QN + CC. Since the NLS part has been treated by Ogawa-Tsutsumi 
[15] and Holmer-Roudenko [8], while the cross terms are higher order, the main 
problem for us is to control QN. Indeed, our way of the localization is motivated 
by a better cancellation in QN, while some other multipliers such as AXv in (13. 2p 
could make the other terms simpler. 

It is further decomposed QN = (QN X + QN 2 + QN 3 )/2 with 

QN! := (v\{XA 1 + A x X)v) = (i^\X(A 1 + A\)u) = (u\Xu), (3.7) 

where we used the symmetry of the bilinear form as well as X = X* and Aq = —Aq. 
Putting 7] := D~ 1 v, the other two terms are computed as follows. 

QN 2 : = (vlDXD-^v) = (r]\D 2 XA^r]) = (Vr]\VXA^T]) 

= (VvlXA^v) + (V V \(Vij R )A. lV ), 
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where we used DA^i = A\D and VA_i = A{V, 

QN 3 := = (ulAiDXD- 1 ^) = {r]\DA l DXrf) = (VrjlA^Xrj) 

where we used DA\D = — V • A{V. Hence 

QN 2 + QN 3 = (V V \X{A 1 + Al)V V ) + (Vjj|(V^)A_h? + A 1 (V4j r )t ] ) 
= (Vr)\XVr)) + 2(Vt 1 \(V^ r )x ■ V77) + (Vr]\r]A d V^ R ) 

= (VrilXVri) + 2(rjr\7i>' R Tj r ) - -(\r)\ 2 \A d - 2 Ai/; R ), 
where we used the radial symmetry of ip R but not of 77. Thus we obtain 
QN = (u\Xu)/2 + (V V \XVri)/2 + (vrW R Vr) ~ \(\v\ 2 \A d ^ R ). 



(3.9) 



(3.10) 



(3.11) 



The first two terms are less than = ||Vr/||| since ip R < 1, while the third term 
is nonpositive since ip' R < 0. The last term is bounded from above and below b}H 



Pr ■-- 



J\x\~R R 2 



dx < ||V?7 



(3.12) 



In short, we have 



QN(t) < \W\\l + 0( PR (t)). 



(3.13) 



Pn{t) ~ >* as R — > 00 for each fixed t, but some uniform decay is needed for the 
main term V^o(t) = 4K + \\v\\ 2 + (1 — c?)(^||m| 2 ) to absorb the error. For that we use 
the equation of 77: 

(id t + aD)r] = D^iidt + aD)(N - \u\ 2 ) = -iD~ l \u\ 2 

and the corresponding integral equation 



t ■ 



(3.14) 



77 = T]° + 77 1 , 77° : = 



iaDt 



—I 



sm(a D(t — s)) 



u(0)| ! 







" /•(*-!)+ 


-./' " 




J(t-1) + 



cos(aD(t — s)) 
5 



(3.15) 



Mi's) | 2 <is 



9 Q 4 

77 + 7/ + ?7 • 



We use the above equation only for very low frequency. More precisely, with a small 
parameter < 6 < 1 independent of t, decompose 77 smoothly in the Fourier space 



(3.16) 



Such an error term does not occur in Merle's virial identity |14j . This is a disadvantage of our 
identity. Nevertheless we can dispose of it using the evolution equation. 



sup PR (t) <||^(0)||^ + ^ 1/2+ [||n(0)||2log(T + 2) + ||z/||l r(0iT;iS) ] 



0<t<T 

+ (fiR) 2 |l Z/ llL t °°(0,T;L2) 



(3.17) 
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then we have ||?/>5||2 < ^"^Mh- For the low frequency part, we have 

\\v <s\\m = \W<s(0)\\^ 

\\vls\\ H ^ + <\\\um 2 \W <\H0)\\l 
\\v<s\\h-^+ ^ a|lhl 2 IUr^ ~ "Il u ( )ll2, 

and by the L°° decay of the wave equation, 

II^IIb-^ < r ^ ^11^)11^ < H0)||>g(t + 1). (3.18) 
Jo \ l ~~ s \ 

Thus we obtain 

WksWl^t^) £ lk(0)||ilog(T + 2) + ||^||i r(0iT;ii) , (3.19) 

and so 



(3.20) 



(3.21) 



Next we estimate the cubic cross terms 
CC = CCi + CC 2 + CC3, 

Cd := -2(uu\A Xu + XA m) = -2(^N4 + XA d )\u\ 2 ), 
CC 2 := ^(u\(XA 1 + A 1 X)\u\ 2 ) = (u\(XA 1 + ri;' R /2)\u\ 2 ), 

CC 3 := -(u\(DXD- 1 A 1 + A 1 DXD~ 1 )\u\ 2 ). 

For the last term we use the commuting relations: 

A X DXD- X = DA^XD- 1 = D(XA_ X + r^ZT 1 

= DXD~ l Ax + Drip' R D-\ 

and so 

CC 3 = (ulXA^ul 2 ) + CC' 3 , 
CC' 3 := (u\([D,X]D- 1 A 1 + [D^' R \D- l /2)\u\ 2 ). 
Hence 

CC = (1 - d){u\\u\ 2 ) + - d)(V> fl - 1) - 3r^/2]|u| 2 ) + CC^, (3.24) 

and the second term on the right is bounded by 

/ \vu 2 \dx< ||f||2||M||2||w|U'»(H>ii) ^ ^"l^lblkll^llVwII^ 2 , (3.25) 

since the functions in the brackets [] vanish on \x\ < R. We used the radial Sobolev 
inequality 

<p{x) = <p(\x\) 6 H\R 3 ) =► |M|^ (R3) < IMlfllV^lf ■ (3.26) 



(3.22) 



(3.23) 
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For the commutator terms CC' 3 , we use the elementary commutator estimate 

\\[D,f]g\\^<mWf)\\ L 4g\\ L2 , (3.27) 
together with the (radial/nonradial) Sobolev 

\\XU 2 \\ 2 < ||ara||oo|M|2 < lkll2 /2 ||Vu||2 /2 , 

WD-'H'h < HM 2 || 6/5 < NI2NI3 < ll^llfllVnllf . 
Since || ^(V^^Hi = CR~ l , we thus obtain 

\CC' 3 \ < WvhR-'iWD-'V ■ x\u\% + WD-'lufM < R-'WvMuff 2 ||V«||f . (3.29) 
In short, we have obtained 

CC = (1 - d)(HM 2 ) + 0(iTiz/|| 2 |M| 2 /2 ||Vw|| 2 /2 ). (3.30) 
Finally we estimate the NLS part 

NS/2 = (-Am - \u\ 2 u\A Xu + XA u) 

= (Vu\\7(A X + XA )u) - (r^JM 4 ) - l$ R \(rd r /2 + d)\u\ A ) (3.31) 
=: iVSi + NS 2 + iV5 3 . 
For the first term NS\ we use 

V(A X + XA ) = A 2 VX + XV A + [V^r]A 

= A 2 XV + A 2 [VV'rJ + XA 2 V + [V^]A (3.32) 
= (A 2 X + XA 2 )V + 2[V^]rS r + [A 2+d V^ R ], 

where the bracket denotes the multiplication with the inside function. Using Aq = 
— Aq as well, we obtain 

NSx = (Vu\2XVu) + 2(u r \i>' R ru r ) + -(V|w| 2 |^ 2 + d V^>. (3.33) 

Since i/) R < 1 and ip' R < 0, the first term is less than 2||Vu|| 2 and the second is 
nonpositive. The last term equals 

-~(\u\ 2 \A d Ai; R ) < \\u\\ 2 2 \\A d AMoo < R^\\u\\ 2 2 . (3.34) 
The quartic terms equal 

NS 2 + NS 3 = ~{{rd r + d)^ R \\u\') = ~\\u\\ 4 4 - ±(pl/ R + d(^ R - l)||n| 4 ), (3.35) 
and the last term is bounded by 

||«|ll4(|x|>H) < lhll2ll«lli-(H>fl) £ ^~ 2 |hll2l|Vn|| 2 , (3.36) 
using the radial Sobolev inequality. In short, we have obtained 

NS/2 < 2K(u) + 0(iT 2 ||u|||||V«|| 2 ). (3.37) 
Gathering the above estimates on QN, CC and NS, we obtain 
V R < 4K(u) + ||i/|| 2 + (1 - d)(u\\u\ 2 ) + 0(p R ) 

+ OiR^WvhWulff 2 1| V«|| 2 /2 ) + 0(ir 2 ||u||||| Vm|| 2 ), 



(3.38) 



Mlif (0,T;L2)] 



(3.39) 
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and 

sup PR <\Wm\l + R- 1/2+ [\\u(0)\\l log(T + 2) + 

0<t<T 

+ 2 \\ u \\l™ (0,T;L|)- 

Also we have 

|VB|<i2[|H| a ||V«|| 2 +||JV||a. (3.40) 
Now suppose for contradiction that 

sup || M (t) + ||iV(t)|| L 2 <Me [l,oo), (341) 

then |M|l2 < M 2 and |Vr| < RM 2 . The variational estimate provides us with an 
upper bound 

Voo = AK(u) + + (1 - d)(v\\u\ 2 ) < — k (3.42) 

for some k ~ J\{Q\) — [E z (v) + X 2 M(u)} > 0. We can first choose < 8 <C 1 so 
small that H^^O)!! 2 . <C k. Secondly we can choose R ^> 1 so large that 

iT 1/3 M 2 log(i?M 2 /(«5K)) < «, (i?" 1 / 3 + (5R)~ 2 )M A < «, (3.43) 

where log(-RM 2 /(<5ft)) may be replaced with (RM 2 /(8k)) 1 / 6 for example. Then for 
< t < RM 2 /8k =: T, we have V R < -k/2, and so \V R (T) -V R (0)\ > kT/2 - RM " 
which is contradicting the above bound on \V R \. 



28 ' 



4. Concentration-compactness procedure 

It remains to prove the scattering in part (a) of Theorem 11.11 Thanks to the 
variational estimates in Section 2, we can proceed as Kenig- Merle. For each < 
a < J(Q) and A > 0, let 

£ x (f,g) := \~ 1 E z (f,g) + XM(f), 

K+(a) := {(f,g) G H l r x L 2 \ £ x (f,g) < a, K(f) > 0}, (4.1) 

S x (a) := sup{||(«,iV)|| s | (u(0),N(0)) e JC+(a), (u, N) sol.}, 

where S denotes a norm containing almost all the Strichartz norms for radial free 
solutions, including Lf(H l x L 2 ). See (I4.25P for the precise definition. For any time 
interval /, we will denote by S(I) the restriction of S onto I. 

From Corollary 12.31 we already know that all solutions starting from )C~^(a) stays 
there globally in time. What we want to prove is the uniform scattering below the 
ground state energy, i.e. S\(a) < oo for all a < J(Q). Let 

E* x := sup{a > | S x (a) < oo}. (4.2) 

The small data scattering in [7] implies that > 0, and the existence of the ground 
state soliton implies that E\ < J(Q). We will prove E x = J(Q) by contradiction, 
and thus finish the proof of Theorem ll.il (a). The main result in this section is 

Lemma 4.1 (Existence of critical element). Suppose E x < J(Q), then there is a 
global solution (u,N) in K x (a) satisfying 

£x(u,N) = E* Xl || (u, N) ||5(-oo,o) = UK A0lko,oc) = oo. (4.3) 
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Moreover, {(u,N)(t) | t G IR} is precompact in H\ x L? x . 

We will prove this lemma by following the concentration-compactness procedure. 
The main difference from NLS is that we need to work with the solutions after 
the normal form transform. In particular, we have some nonlinear terms without 
time integration (or the Duhamel form). Besides that, we have various different 
interactions, for which we need to use different norms or exponents. 



4.1. Profiles for the radial Zakharov. First we recall the free profile decompo- 
sition of Bahouri-Gerard type [TJ. Actually we do not need its full power, as we can 
freeze scaling and space positions of the profiles thanks to the radial symmetry and 
the regularity room of our problem. Hence the setting is essentially the same as the 
NLS case 0. 

Lemma 4.2. For any bounded sequence (f n ,gn) in H} x L? r , there is a subsequence 
(fni9n)> J £ NU {oo} ; a bounded sequence {P , g J }i<j<j in H} x L? r , and sequences 
{tn}neN,i<j<J ^ such that the following holds. For any < j < J < J , let 

-itA fl at f-i\ ._ JtaD l 

Mi 



u n (t) := e-^C N n (t) := e ltalJ g' n 

ui(t) : = e - J (*-*") A F, N{(t) := ^"^"V, 
J J 

N> J :=N n -^l 

j=l 3=1 



u> 3 : = 



(4.4) 



Then for any j, k 6 {1 . . . J}, we have t^ := lim^oo t 3 n G {0, ±oo}, 

j ^ k lim \ti -t k n \ = oo, (4.5) 



(4.6) 



(w> J , iV^ J )(4) -> weakly in H 1 x L? as n ->■ oo , 
(u> J , iV> J )(0) ->■ wea% m if 1 x L 2 as n ->■ oo, 

and 

lim limsup[||u> J || 1/2-^ + |iV n >J | B - 3/ 2- S R - 3 /2 +i J = 0. (4 7) 

Remark 4.1. 1) fl4.5p - fl4.6p implies the linear orthogonality 



lim \\u n (0)\\m - lk(0)llln - 11^(0)11^ = 0, 

J 

lim M( Mn (0)) - ]T MK(0)) - M( M > J (0)) = 0, (4.8) 
lim IK(0)||| 2 - £ ||Ni(0)||| 2 - ||iV> J (0)||i 2 = 0, 
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as well as the nonlinear orthogonality 

lim K(0)|| 4 L4 " £ lk(0)||t. - lk> J (0)|| 4 L4 = 0, 



lim BsWO)) - E^K(O)) - £ 5 (u> J (0)) = 0, 

i=i 
J 

lim KK(o)) - E^K(o)) - k(u>\o)) = o, 

3=1 

J 

lim E z (u n (0),N n (0)) - £^(<(0),N£(0)) - £ Z « J (0), iV> J (0)) = 0. 



(4.9) 



3=1 

The same orthogonality holds also along t = t J n instead of t = 0. 

2) The norms in (14. 7p are related to the Sobolev embedding L? C -B^ 2 . Inter- 
polation with the Strichartz estimate extends the smallness to any Strichartz norms 
as far as the exponents are not sharp either in LP or in regularity (including the low 
frequency of N). 

We call such a sequence of free solutions {(u J n , N{)} nS N a free concentrating wave. 
Now we introduce the nonlinear profile associated to a free concentrating wave 

(u n (t),N n (t)) = U(t-t n )(f,s), too = lim t n G {0, ±oo}, (4.10) 

where U(t) = e~ ltA © e ltaD denotes the free propagator. With it, we associate the 
nonlinear profile (u, 9T), defined as the solution of the Zakharov system satisfying 

(u,N) = U(t)(f,g)+ [ U(t- s)(nu,aD\u\ 2 )(s)d Sj (4.11) 

which is obtained by solving the initial data problem (if t^ = 0) or by solving the 
final data problem (if = ±00). When too = ±00, the existence of wave operators 
will be given at the end of this paper as appendix . 

We call (u n (t), %ln(t)) '■= (u(£ — t n ),D f t(t — t n )) the nonlinear concentrating wave 
associated with (u n (t), N n (t)). By the above construction we have 

||K,N n )(0)-(u n ,9T n )(0)||^ xL2 

= \\( u ,m)(-t n )-u(-t n )(f,g)\\ mxL 2^o. [ - ] 

Given a sequence of solutions to the Zakharov system with bounded initial data, we 
can apply the free profile decomposition Lemma 14.21 to the sequence of initial data, 
and associate a nonlinear profile with each free concentrating wave. If all nonlinear 
profiles are scattering and the remainder is small enough, then we can conclude that 
the original sequence of nonlinear solutions is also scattering with a global Strichartz 
bound. More precisely, we have 

Lemma 4.3. For each free concentrating wave (u{,N{) in Lemma\4^ let (u J n , 9T{) 
be the associated nonlinear concentrating wave. Let (u n , N n ) be the sequence of 
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nonlinear solutions with (u n ,N n )(0) = (/„,(?«)■ If ||( u ra>^n)lls(o,oo) < 00 for all 
j < J, then 

limsup || (u n ,N n ) || 5(0,00) < 00. (4.13) 

To prove Lemma 14. 3\ we need some global stability. In the next subsection, we 
will refine the normal form reduction and the nonlinear estimates that was used in 
[7J , and then prove Lemma 14.31 and Lemma 14.11 

4.2. Nonlinear estimates with small non-sharp norms. In order to obtain the 
nonlinear profile decomposition, we need that the non-sharp smallness (14 ,7p is suffi- 
cient to reduce the nonlinear interactions globally. The idea is to use interpolation, 
thus we need to do some refined estimates than in [7J , more precisely, to avoid using 
the sharp (or endpoint) norms with L\ or L£°. 

4.2.1. Modifying the nonresonant 'part. The first problem in following the Strichartz 
analysis in [7J is the L^-type norms. In fact, one can observe that the use of L^-type 
Strichartz norm for N is inevitable for the low-high interactions of nu in very low 
frequencies, since the regularity exponent becomes bigger than that for the dual 
Schrodinger admissible exponent as we move the Strichartz norm of N to L 2+ . 

However, this problem can be avoided by applying the normal form to those 
interactions. In fact, there is no resonance in very low frequencies because 

-|£| 2 ±a|£-r/| + M 2 ~a|£-?/| (4.14) 

when all of |£|, |£ — r)\, \r)\ are small. Hence we include them into the "non-resonant" 
interactions, which are integrated in time before the Strichartz estimate. 

The second problem is that our solution is no longer small, so the nonlinear terms 
without time integration (i.e. the boundary terms from the partial integration) 
do not contain any small factor for the perturbation argument. To overcome this 
difficulty, we shrink the "non-resonant" part to either higher or lower frequencies, 
for which we gain a small factor, depending on the frequencies, from the regularity 
room. Hence our decomposition into the "resonant" and "non-resonant" interactions 
depends on the solution size. 

Thus we are lead to divide the bilinear interactions nu and |w| 2 as follows. Let 
u = J2k& ^ u ^ e the standard homogeneous Littlewood-Paley decomposition such 
that supp FPkU C {2 fc ~ 1 < |£| < 2 fe+1 }. For a parameter f3 > 5 + | log 2 «|, let 



XL 


:={(j,k)eZ 2 


j > max(fc + 5, /3)}, 


RL : 


= {(j,k)ez 2 \ 


\j\< P and j > k + 5}, 


LL : 


= {(j,k)ez 2 \ 


max(j,k) <-/?}, 


LH : 


= {(j,k)ez 2 


k > mm(j - 5, -/?)}, 


HH 


:={( 3 ,k)eZ 2 


\j — k\ < 5 and max(j, k) > f3}, 


RR : 


= {( 3l k)ez 2 


max(j,fc) < /?}, 



(4.15) 



and LX := {(k,j) | (j, k) e XL}. Then 



Z 2 = (XL U LL) U (RL U LH) = (XL U LX) U (HH U RR), (4.16) 
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where all the unions are disjoint. For any set A C Z 2 , and any functions f(x),g(x), 
we denote the bilinear frequency cut-off to A by 

(fg) A = F- 1 f VaKC ~ V)9(v)dv := E ( P M P <*9)- (4.17) 

(j,fc)GA 

For the nonlinear term nu, we apply the time integration by parts on XL U LL, 
where the phase factor oj = — |£| 2 ± a|£ — 77 1 + |?7| 2 is estimated 

M~ie-^i(e-^)~ie-^i(e>, (4.18) 

which is gained in the bilinear operator 

n± (/, S ) := I P XLULL /« 7 ^ 

7 -Kr±«K-^l + l^r (4.19) 

fi(/,s):=^{fi + (/,s) + fi_(7,s)} 

For the nonlinear term mw, we integrate by parts on XL U LX. Then we get a 
bilinear operator of the form 

7 If - w 2 - M 2 - 

After this modification of the normal form, we can rewrite the integral equation 
for (11.41) as follows. Let 

u := (u, N), v° := U(t)u(0) = (e" itA M(0), e itaD N(0)). (4.21) 

For the fixed free solution u°, the iteration v! y u is given by 

u = if - U(t)B(u(0), u(0)) + B(u, uf) + Q(u', u') + T(u', &, &), (4.22) 

where the bilinear forms B, Q and the trilinear form T are defined by 

B{u x ,u 2 ) :={Sl{NuU2),DVt{y, u U2)), 

Q(Ul,U 2 ) ■= / U(t - s)((n 1 U 2 ) LHUaL ,D(u 1 U 2 ) HHUaLuLa )(s)ds, 

Jo 

T(ui,U2,Uzl) := / U(t - s)^^^^),^) + Q(N 1 ,n 2 u 3 ),Dfl(u 1 ,n 2 u 3 ))(s)ds. 



For brevity, we denote 

NL(u u u 2l u 3 ) :=B(u l ,u 2 ) +Q(u 1 ,u 2 ) + T(ui,u 2 ,u 3 ), NL{u) := NL(u,u), 

B{u) :=B(u,u), Q(u) := Q(u,u), T{u) := T(u,u,u). 

We can estimate each term in the Duhamel formula using some powers of Strichartz 
norms with non-sharp exponents. For brevity of Holder-type estimates, we denote 
the space-time norms by 

(b,d,s) :=L] /b B[ /d ^ 

(6, d±e, s)+ := (b, d + e,s) + (b,d- £, s), (4.23) 
(b, d ± s, s) n := (b, d + e, s) H (b,d — e, s). 
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Using the above notation, we introduce nearly full sets of the radial Strichartz norms 
for the Schrodinger and the wave equations. Fix small numbers 

< k < e < 1, (4.24) 

and let 

SS:= (D)- l [(0,l 



o) n 


(11- 


K 


2 


v 2' 10 


3' 


5 ~ 


l l 


K 1 




k), 


2'4 


3' 4 






by using 








1 \ n \\ 


3 


-« 


Lf (B, 


X) 


+B, 



sw := (o,^,o)r 



(4.25) 



M Ly := M ; \\ n \\Y '■= \\ n \\ - 3 -«5 - 3 +6 i Z := I X 7. OfA 

In the nonlinear terms, we should choose appropriate Strichartz exponents so that 
all can be controled by interpolation between S and Z. For that purpose, we will 
choose (b, d, s) for u and N respectively to be H s admissible with < s < 1 and 
L 2 admissible for radial functions. Moreover, b < 1/2 and (b,d) ^ (0, 1/2). Besides 
that, we will use the sum spaced with small e > for N and the intersection for 
u, so that we can dispose of very low or high frequencies, and sum over the dyadic 
decomposition without any difficulty 

4.2.2. Bare bilinear terms. First consider the bilinear terms which do not contain 
the time integration, namely the boundary term in the transform. In the equation 
for u, Q(n,u) is roughly like (D)~ l {D~ l n)u for each dyadic piece. 

Lemma 4.4. (a) There exists 9 > such that for any N and u, we have 

Mn,u)\\ L ^ <2- e >||y||n||y||n||^|M|^ (4.27) 

||n(n,«)||55<2- tf/, ||«||^l|nKl|n||y|l«lli- (4-28) 
(b) There exists 9 > such that for any u and v! , we have 

II^K^OIlL^^lkliyikll^lkll^ll^lll, (4.29) 

||^KnO||w<2^||n||y||n'||y||n||^||n'||^. (4.30) 

Proof, (a) By the Coif man- Meyer-type bilinear estimate on dyadic pieces (see 
Lemma 3.5]), we have for (j, k) G XL, 

llnfaj.UfcJHLooffi < ||-D~ 1 nj|| ( o,i ±e)0 ) + ||Wfc||(oi± e ,o)n 

< O -0/10 1 1 r)-l„ II . II.. II . 

and for (j, k) G LL, 

\\^{nj,U k )\\ L oo H i < H^ -1 ^!! (0i ^ ±£j o) + ll^fcll (0,|l±e,0)n 

< 2~ /?/1 ° WD^Uj || (o, i ± e ,l) + IKH (0,ll±£,0) n • 

Since the right hand side is bounded by ||n||£,oo£2||it||£ooffi via non-sharp Sobolev 
embedding, we obtain, after summation over dyadic decomposition, 

||^(n^)|U^ 1 <2^/ 10 || U ||y||n||y||n||^||«||^, (4.31) 



4 This is because N{0) G L 2 while u(0) e H 1 = L 2 n H 1 . 
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for some small 9 > 0. Similarly we have, for (j, k) e XL, 

\\^{nj,u k )\\ {D) -i { i t 3__, t 2_ K) < \\D- 1 n j \\ { i^7__. ±£ ^ K)+ \\u k \\ { i^ ±£)0)n 

<2-^lD-S.|| (i x_ f±£ ,^ )+ ||^|| (| ^ ±£ , 0)n , 
and for (J, k) G LL, 

||^K-,« fc )|| (D) -i ( i^_ f)0) < || J D~ 1 n,-|| ( i ) ^_ f±£)0)+ ||« fc || ( i ) ^ ±£!0)n 

~ II J ll( IlT 5— 3±e ,-2q)+ II «ll(4,3o±E,0)n 

Hence in either case we can control by non-sharp norms, so 

\\n(n,u)\\ ss < 2-^ 20 || W ||fe1n||y \\n\\ d Y \\u\\ d x . 

(b) We may assume (j, k) e XL, since the other case LX is treated in the same 
way Similarly to the above, we have DQ(fj,g k ) ~ 1 (fj9k)i so 



(4.32) 



(4.33) 



(4.34) 



||-D^(WjX)IU°°£ 2 ~ \\( D ) D ^(Uj,U k )\\ L ^ {L 2 +L 6/ 5) 



< 



Ml(0,i± e ,0) + IKJ' 



(0 i± £ ,0)+ 



<2-^°||« 



j II (0 ±±e i)+ II "fell (0,i±E,0)+ > 



hence 



Similarly, 



\\DQ(u,u')\\ L ^ <2-^°\\u\\]r s e \\u% 



ss 



\u\ 



X 



u 



11^(^,011(1 i » i „, < \\(D)Dn( Uj ,u' k )\\a l0) 



< 



2-^ 10 \\ Uj \ 



(i,ii)||%||(i,i,o), 



and so 



\\DQ(u,u')\\ sw 
Thus the proof is completed. 



<2-^ 10 H«||^||« / || 1 c 



\ss 



u\ 



X 



II 



x ■ 



(4.35) 

(4.36) 

(4.37) 

(4.38) 
□ 



4.2.3. Duhamel bilinear terms. Next we consider the remaining bilinear terms in 
the Duhamel form after the normal form transform. Here we have to use the radial 
improvement of the Strichartz norms. For brevity, we denote the integrals in the 
Duhamel formula by 



/«/:= t e- t{t - s)A f(s)ds, I N f:= f e ^ aD f(s)ds. 
Jo Jo 



(4.39) 



Lemma 4.5. (a) There exists 9 > and C(/3) > 1 such that for any N and u, we 
have 



\\h{nu) LH \\ ss <C(m^\Hlw\H\Y 



\u\ 



X ! 



\\iu{nu) RL \\ss <c{mA\~s e \H\t\HY\Hx- 

(b) There exists 9 > and C(/3) > 1 such that for any u and u' , we have 

\\InD{uu') hh \\sw <cmHts\\Ats\Hx\\Ax, 



\\i n d{uu') rr \\sw <cmn l ssh 



■'^WvW 9 \\v'\\ e 
\ss II "llxll" \\x- 
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Proof. In this proof we ignore the dependence of the constants on /3. 
(a) For (j, k) £ LH , we have for < s < 1, 

ll n i M fc|l(l-2 £ ,i+2 £ , S +2 e ) < |kj||(i_ £! | ± |-i_ £ ) ± |kfc||(i_ e) i + 2 e ±f, s+ | +3e )n , 44Q . 

K-|l(±- e ,i±§,-±- £ )J M fc|l(i- £ ,i+2 £ ±§,f+3 £ )n' 



where in the second inequality we used that k is bounded from below. Since the 
left hand side is if s -admissible norm for the Strichartz estimate (without the radial 
symmetry), we obtain the full Strichartz bound in H 1 . 

For (j, k) £ RL, we may neglect the regularity of rij and the product, since 
their frequencies are bounded from above and below. Using the radial improved 
Strichartz, the full H 1 Strichartz norm is bounded by 

K-U fc || { ! +2ei 3_ 3ei0) < IMI(I_ e ,i,0)IKII(3 £ ,|-3 £ ,0)- (4.41) 

Summing these estimates over dyadic pieces in the specified regions, and using 
non-sharp Sobolev embedding and interpolation, we obtain 

\\Iu(nu) LH \\ ss < IMIs5 e IMIwH™llylMlx, 



\Iu(nu) RL \\ ss < \\u\\l 



U- 



SS \\ n \\sw\\ n \\Y \\ u \\x- 



(4.42) 



(b) We consider only the case j > k for Uju' k , since the other case is treated in 
the same way. For (j,k) £ HH, 

||%<|| (1 _ e> i +fejl+e) < ||^|| ( i_ f: i +f| i +f) ||<|| ( i_ f) i +fj i +f) , (4.43) 

and in the case (j, k) £ RR, since j is bounded from above, 

IK' u fcll(!+e,f ,f+e) ~ ll u j'll(i-e ,|+2e,|+e) ll M fc II (2 £ i-2e,0) 



Hence 



w2 1 '4'4 1 ' v 2 '4 1 '4 1 ' " v ' 2 

~ ll M ill(|-e ,\+2e ±)\\ U k\\{2e ,|-2e,0)- 



i-D(W) ra ii(i- £ ,i + | £)£ ) < ikiiyh'iiyihiixik'iix, 

'Hi-flu . Me n 'i 



|D(W) QjL+Lq || ( i +£! | _i +£) < llullg/l 



(4.44) 



(4.45) 



M 1155 ll' u llxll M NX- 



The left hand sides are L 2 -admissible norms for radial functions. Thus the proof is 
completed by the radial improved Strichartz. □ 

4.2.4. Duhamel trilinear terms. Finally we estimate the trilinear terms which appear 
after the normal transform. These are supposedly the easiest, but there is a small 
complication due to the fact that we have to use negative Sobolev spaces for N in 
some of the products: 

ra*3 s ii/ii^iisHbj, 

< s < 3/q, 1/r = l/p + 1/q- s/3. 
In the next lemma, the constant may decay as /3 — > oo, but we do not need it. 
Lemma 4.6. (a) There exists 9 > such that for any u,v,w,n,n', we have 



\I u Q(D(uv),w)\\ss <\\u\\\- 6 

\sw 



ss \\ v \\ss \\ w \\ss iFllxlPllxlFllx- 
|/„n(n,n / «)|| S5 <||n||^||n / ||y||u||y||n|ltlln / lltllur fl 
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(b) There exists 9 > such that for any n,u,u', we have 
\\I N DQ(nu,u')\\ sw +\\I N DQ(u,nu')\\sw < ||n|lswll u llssl u 1l5s e ll w llyll u llxll u 'llx- 
Proof, (a) Since Q(D(uv)j,w k ) ~ (D)~ 1 ((uv)jW k ), 

\\n(D(uv)j,w k )\\ L i H i < ||w|U3L6 11^11^6 ||w|| L 3 X 6, (4.47) 
and by non-sharp Sobolev embedding and interpolation, 

\\n(D(uv),w)\\ L i m < Ik||^lkll5^lhliyil«llxlkllxlhllx- (4-48) 
For Q(rij, (n'u)k), we have either 2 J ^> 2 fe or 2 J + 2 fc < 1. In the first case, we have 
\M n] ,(n'u) k )\\ Lim < \\D- 1+Se nj \\ { i_ e ,2 e ± § ,o) ± P~ 5e («'")fcll(i+ e ,|-2 e ±|,o)n „ , , 

< II II II 'II II II ( 4 - 49 ) 

~ ll n ill(|-e,2e±f ,-l+5£)±ll ?7 ' II (2e,|-|e±f ,-5e)± II M ll (|-e,2s±f ,5e) n ' 

where we used the product estimate for negative Sobolev spaces for n'u. In the 
second case 2 j + 2 k < 1, we have 

\\n{nj, {n'u) k )\\ L i H i < ||-D~ 1 n i ll(J- e ,f±|,o) ± ll( ra ' u )*ll(^+ e ,J-§±§ 



< ||nj|| ( i_ £) | ± | _!) ± ||n || (2£) i_7 e± | _ 5£ ) ± ||w|| ( i_ £) u e± | j5e ) n . 
Hence, by non-sharp Sobolev embedding and interpolation, 



(4.50) 



\\n(n,n'u)\\ Lim < llnliyiln'liyilwIiyilnll^lln'II^IKII^. (4.51) 
(b) We have DQ ~ (D)" 1 on each dyadic piece, so 

\\Dn{{nu)j,u' k )\\ L i L ^ < ||fi((nu) i ,M , fc )|| (1) 5_| e _ 5e) 

£ ||Mi|| ( i +£) |_ £ _ 5£ )lkfell(i- £ ,i+ £ ,5 £ ) (4-52) 

~ lkll(2£,i-|e±f _5e) ± ll u ll(i-e,I+e±|,5£) n ll W fell(|-e,|+e,5e)' 

where we used the product estimate twice, but did not use any restriction on j, k. 
Hence we have the same estimate on Q(uj, (n'u) k ), and so 

\\Dh(nu,u')\\ L i L 2 + \\DQ(u,nu')\\ L i L2 

^IHI^IMI^IMI^IHI^IMI^IMIi- 

Thus, the proof is completed. □ 

Note that in the above estimates we needed the L£°-type norms only for the bare 
bilinear terms, but not for the Duhamel terms. Thus we have obtained 

Lemma 4.7. There exist 9 > 0, rj > and C(/3) > 1 such that for each (3 3> 1 and 
any ui,U2,u^, we have 

2^\\B(u u u 2 )\\ s + \\Q(u u u 2 )\\s/C(P) < !|«i|ini«2|ir e H«illill^|||, 

wn^u^uMs < ii^iirii^iirii^iirii^iilii^iizii^iii. 

For the Duhamel terms we have also 

\\Q@i,b)\\s<cm\M§\\Ms> (4 55) 

||r(ui,ii2,M3)||s < ll«l|lsll«2|ls||M3|l5, 
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where 

S :=SS x SW, 

^■^{Dr^l-y^m 1 -,^-! 2 --.)], ( 4.56) 

SW := (77, 77, --77) n (-, , k). 

w ' 2 2 ' 4 /; T 4 3 4 ; 

4.3. Nonlinear profile approximation. We will prove Lemma 14.31 by the follow- 
ing two lemmas. 

Lemma 4.8 (Stability). For any A > and a > 0, there exists q > with the fol- 
lowing property: Suppose that u a satisfies \\ u a \\ 5(0,00) — A an d approximately solves 
the Zakharov system in the sense that 

u a = U(t)u a (0) - U(t)B(u a (0)) + NL(u a ) + e 

and ||e||s(o,oo) < Then for any initail datau(0) satisfying ||u(0) — u a {fS)\\ H i xL 2 < q, 
there is a unique global solution u satisfying \\u — u a ||s(o,oo) < °~- 

Proof. Denote = u a — u, then ||"U>(0)||#i X £2 < q and 

=U{t)u^) - U(t)B(u a (0)) + NL{u a ) + e + U(t)B(u(0)) - NL{u). (4.57) 

Thus 

||«>||s + \\B{u a ) - B(u)\\ s + \\Q(u a ) - Q(u)\\ s + \\T(u a ) - T(u)\\ s . (4.58) 
Noting that Z D S, by (14. 54ft we have 

\\B(u a ) -B(u)\\ s < \\B(u a ,u > )\\ s + \\B(u > ,u a )\\ s + ||5(m > ,t2 > ) || s 



< 2^A\\u > \\ s + 2- 



6P\\rf ||2 



■a 



>\\s- 



By f)4.55p . we have 

\\Q(u a ) - Q(u)\\ s <C{p) {\\uaWsWM\s + Kill) 

I ~ 
Is 



\\T{u a ) - T(u)\\ s <C (||« |||||«>||s + IKIIslKlls + Kill) 



So 

Klls < 2? + (2- e ?CA + C(/3)\\u a \\s + CIKIIDKIU 

+ (2~^C + C(/3) + CII^II^II^IH + CK|||. 



(4.59) 



(4.60) 



Choose /3 = (3(A) such that 2 e/3 CA < j. Then we subdivide the time interval 
[0, 00) into finite subintervals Ij = [tj, j = 1, ■ ■ • , J, J = J (A, (3) such that 

^09)11^.115(1,) + ^IKII|(/ J) <^ ( 4 - 61 ) 

for each j. Let q = q(A, a, 0, J) small such that 

C(/3)8 2 ^<1, 8 2J q<^a. (4.62) 
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Then by ( 14.60P on Ji, we have 1 1 w t> 1 1 ^(Ti ) < 8? an d 

||«>(^2)||h1xL 2 

<\\U(t 2 - ^u^Wh^ + ||tf(i 2 - t^Biuait^)) - U(t 2 - t^Biuit^Wn,^ 

+ \\B(u a (t 2 )) - B{u{t 2 ))\\ mxL 2 + \\Q(u a ) - Q(u)\\ sih) 

+ \\T(u a )-T(u)\\ S{h) + \\e\\ S{Il) 
<2q + 4 • 8? < 8\. 

Using the same analysis as above, we can get ||w>||sf(7 2 ) < 8 3 ^. Iterating this for 
I 2 , 7 3 , . . . , Ij, we obtain ||« a — ui\\s ^ 8 2J ^ ^ 17 , the desired result was obtained. □ 

With J close to J and large n, our approximate solution is given by 

j 

< = « N$ := ]T«, mi) + (u> J , N> J ). (4.63) 
i=i 

To prove Lemma [4.31 we only need to prove that u n is an approximate solution of 
the Zakharov system. In fact, we have 

Lemma 4.9. Suppose that ||(u£, 9l£)||s < oo for all j < J, then 



j 

lim hmsup \\U(t)B(u J n (0)) - NL(u J n ) - £>(t)B(0* (0)) - iVL(<)]|| s = 0. 

J— >J n— )-oo . -, 



Note that ||(u^, O^Hs does not depend on n. 
Proof. By triangle inequality, it suffices to show that 
hm || 5>(*)2J(t&(0)) - NL(ui)} 



n— >-oo 



(4.64) 



(4.65) 



- u£(0)) - nlC£K)\\\s = o, 

and 

lim limsupH [[/(*)£«(())) - NL{u J n )} 

n— >oo 

- [tf(f)S(t#(0) - u> J (0)) - iVL« - u> J )]\\ s = 0. 
In fact, 

L.iT.5o/gSl<^(||i?(<,<)|| 5 + ||Q(<,<)|| 5 )+ £ ||T«,<,^)|| 5 . 

i^j i^j or jjtk 

For each i ^ j, we have |^ — ^| — > oo. for the subsequence f n — P n — > oo, we have 
by 

I|£(<,#JIU < ||B(u*(- -<),#(• -ii))IU ( _ 00j(4+ ^ )/2)nS((4+ 4 )/2i00) 

r i ( 4 -66) 

< lliTll^llfP'll 1 - 9 lliPII 9 4- WWW 6 

~ ll U IIS H U IIS [H U Hz(-oo,4-t?J/2)ll U + l|U ^ l|U llz((f n -^)/2,oo). " 

For each j, by the scattering of u J , 

T lim o l|u- ? |U(|t|>T) = 0, (4.67) 
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so from the above estimate 

||fi(<,<)|| 5 ->0, (4.68) 

as tf n — t J n — > oo. The case f n — P n — >• — oo is treated similarly, as well as the other 
terms Q and T. Thus we obtain 

\\B(&,%)\\ 8 ^0fori?j, 

\\Q(x? n ,ui)\\ s ^0 for i^j, (4.69) 

||T(i£,u*,u£)lls -> for i ^ j or i = j ^ k, 

from which (14.641) follows immediately. 

In order to prove (14.651) . we need a uniform bound on the approximate solutions 
u J n for J J. Note that ( 14. 9 p implies that || (0) x z, 2 1 except for a bounded 
number of j. Let A be the set of j in the latter case. Then for all j (jL A, the small 
data scattering implies that 

\\(<,<K)\\s< ||(ui(0) l Ni(0))|| fllxZS «1. (4.70) 
Then by the orthogonality in H\. x L? x and \ f n — P n \ — > oo, we deduce 

II E^Hs S E ll^lll ~ E Kllfllx^ £ !■ (4.71) 

i^A j£A j&A 

Since the number of the remaining components j & A are bounded, we obtain 

supsup \\u J n \\ s < oo. (±72) 

j<j « 



The left hand side of (14.651) is bounded by 



\B{u J n ) - B{u J n - u>J)\\ s + \\Q{ui) - Q{u J n - u>')\\ s 



(4.73) 



+ \\T(<)-T(u J n -u> J )\\ s . 
By gSD, gZD and ([122]), 

lim_limsup ||5«) - B« - w> J )|| 5 = 0. (4.74) 

One can estimates Q and T similarly. Then ( 14.651) was proved. □ 
Proof of Lemma \4.3[ By the construction of u J n , 

lim ||<(0) - u n (0)\\ H i xL 2 = 0. (4.75) 

By Lemma 14.91 and Lemma 14. passing to a subsequence if necessary, we obtain 
||«n — UnWs ^ 1 fo r large J and n. □ 

Proof of Lemma \4^1 By the definition of E x , there is a sequence of global solutions 
(u n ,N n ) in /Cj~(a) such that 



lim £ x (u n ,N n ) = E* x , lim ||(u n ,iV n )||s(-oo,o) = hm ||(it„, A^ n )||s(o,oo) = 00. 

n— >oo n— >oo n— >oo 



To see this, first note that since (u n , N n ) are in )C x (a), they are bounded in H 1 x L? 
by the energy. Hence if || (u n , N n ) \\$ —> 00 then the L? t part must diverge, and we 
can translate (u n , N n ) in t so that the norm diverges both on (—00, 0) and on (0, 00). 
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For the sequence (u n (0), N n (0)), we use the linear profile decomposition. For the 
associated nonlinear profile (u^, we must have i^(u{(0)) > for each j. In fact, 
if we denote 

i„ n i\ ,._ , l2 , a 2 ll2 



G A fo>) = JaM - ^ = ( ^ - g ) II Vm H2 + y NI2 > 0, (4.77) 
then 

J\(Q\) = inf{ JaM | <^0, tffa) = 0} 

= inf{G A M I <f + 0, tf(^) = 0} (4.78) 
= inf{G A (^) I tp + 0, tffa) < 0}. 

By the orthogonality, 

En" GaK(0)) = lim" f Vg a (<(0)) + G a « j (0)) ] < \E\ < J X (Q X ). (4.79) 

\j"=l / 

Hence, for n sufficiently large, G\(u 3 n (0)) < J\(Q\); and then by the third line of 
< K7m . K(ui(0)) > 0. Noting that 

j 

lim AK(0),iV n (0)) - A«(0), W n (0)) - A« J (0),iV n >J (0)) = 0, (4.80) 

3=1 

we have 

j 

V A(<(0), W n (0)) < hm AK, iVn) = £5. (4.81) 

' n— >oo 

If A(<(0), W n (0)) < #J for all j < J, then we have ||(u£, 9^)|| s < 00 for all j, 
and so by Lemma 14. 3[ 

limsup ||(it„,iV n )|| s < 00, (4 82) 

which contradicts lim^oo \\(u n , A^ n )|| S ( 0iOO ) = 00. Thus, we must have one j < J 
such that 

A«(0),^(0)) = ^. (4.83) 

Without losing generality, we may assume j = 1. Comparing this with ( 14.8ip . we 
have 

K(0),iV n (0)) = Ui-QifW) + (^(O^A^O)) (4.84) 

and 

||« 1 (0),iV> 1 (0))|| Hlxi2 ^(^(O),]^)) ->0. (4.85) 
If t n — > —00, then we have 

llf/^-^X^g^lkcoo^O, (4.86) 

and hence 

||C/(0K(0),iVn(0))|U(0,oo) 

< ||t/(t-t n )(f 1 ,g 1 )|U ( o, oo) + ||( M > 1 (0),iV> 1 (0))b 1 xL 2 ^0. 
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By Lemma [4.71 

\\U(t)B((u n (0),N n (0)))\\ si o,oo) + \\NL(U(t)(u n (0),N n (0)))\\ sio>oo) -> 0. (4.88) 

Then using Lemma l4~8l (with u a := U(t)(u n (0), N n (0)) and (w n (0), N n (0)) as the 
initial data), we obtain 

lim || (u n ,N n ) || s( 0) oo) < oo. (4.89) 

which contradicts || (u n , N n ) \\s(o,oo) — > oo. 

If t n — > +oo, the argument is similar and we obtain a contradiction by using 

||(ttn,^Vn)||s(-oo,0) -> OO. 

So, the only case left is t n — > 0. In this case, 

||(n n (0),iV n (0)) - (f 1 ^ 1 )^^ 0. (4.90) 

Let (it, iV) be the global solution with initial data (w(0), N(0)) = (f^g 1 ), then 
$\{u,N) < E x . By stability, we must have 

||KA0lk-oo,0) = || (U,N) || 5(0,00) = oo, (4.91) 

since otherwise (u n ,N n ) should be bounded either in S(— oo,0) or in 5(0, oo). By 
the definition of E* x , £ x {u, N) > E* x and hence <f A (u, N) = E* x . 
Since (u, N) is locally in S, for any t n G R, we have 

||KA0|| 5 (-oc, t „) = oo= ||KiV)||s(t„,oo)- (4.92) 

Applying the above argument to the sequence (u n (t), N n (t)) := (u(t + t n ), N(t + t n )), 
we see that (n(t + t n ),N(t + t n )) is precompact in H 1 x L 2 . Thus we obtain the 
desired result. □ 

5. Rigidity Theorem 

The main purpose of this section is to disprove the existence of critical element 
that was constructed in the previous section under the assumption E x < J(Q). The 
main tool is the spatial localization of the virial identity. We prove 

Theorem 5.1 (Rigidity Theorem). Let (u, N) be a global solution to (ll.4j) satisfying 
K(u) > 0, and Ez{u,N) + X 2 M(u) < J\{Q\) for some A > 0. Moreover, assume 
{(u, N)(t) : t G R} is precompact in H 1 x L 2 . Then u = N = 0. 

Proof. By contradiction, we assume (it, iV) ^ (0,0). Then by the compactness we 
may assume further u ^ 0, since otherwise iV would be a free wave and dispersive. 
We divide the proof into the following three steps: 

Step 1: Energy trapping. 

We claim that 

c:=MK(u)>0. (5.!) 

If not, then there exists {t n } with t n — > t* G [—00,00], and K{u(t n )) — > 0. By 
the precompactness of {u{t) : t G R}, we get that up to a sequence (u(t n ), N(t n )) 
converges to some (f,g) in H 1 x L? . Then we have K(f) = 0, J\(f) < Ez(f,g) + 
X 2 M(f) = E z {u, N) + X 2 M(u) < J\{Q\). By the variational characterization of Q\, 
we get / = which contradicts to the M(f) = M{u) ^ 0. 
Step 2: Uniform small tails. 
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Let v = dtN — \u\ 2 = n — \u\ 2 . We claim that for any e > 0, there exists R > 
such that at any t G R, we have 



J\x\>R 



(|Vn| 2 + \u\ 2 + m 4 + |u| 6 + \u\ 2 + \d~ x Vv\ 2 + K^r-)dx < 



Indeed, since {(u,N)(t) : t G M} is precompact in iJ 1 x L 2 , by Sobolev embedding 
and the L p -boundedness of -D _1 V, we get that {u(t)} is precompact in L 2 ,L 4 ,L 6 , 
{D~ 1 N(t)} is precompact in if 1 , and {D~ 1 VA r (t), z/(t), D _1 Vz/} is precompact in 
L 2 . Then the claim follows immediately. 

Step 3: Contradiction to the local virial estimates. 

We recall the local virial estimates obtained in Section |3j For any R > 

V R (t) := (Ju\(AX + XA)u). (5.2) 

where X = X* be the operator of smooth trancation to \x\ < Rbj multiplication 
with ipji(x). From the proof in Section [3] and Corollary 12.31 we have 

\V R (t)\<R[\\u\\ 2 \\Vu\\ 2 +\\N\\ 2 ]<R. (5.3) 

On the other hand, from Step 2, Step 1, and Lemma [2.41 we get 

V R (t) =(u\Xu)/2 + (Vr)\XVr))/2 + (r/ r |r^ r ) - -I^A^r) 



1 



2{u\\u\ 2 ) +0( j R~ 1 ||z/|| 2 ||m||2 / ||Vm 



3/2|,„ 1 1 l/2x 



2 



+ (Vn|4XVn) + 4(u r \tp' R ru r ) 

- 3\\u\\l + 0( J R~ 2 (||n|| 2 + ||«||i||Vu|| 2 )) (obtained in Section[3]) 

=AK(u) + \\u\\ 2 2 - 2(^||n| 2 ) + o(l), R ->• oo (by Step 2) 
2 

>(1 i=)K(u) + o(l) > c/2, i? > 1. (by Step 1 and Lemma E2D 

v 6 

Thus we get 

VrW > Vr(0) + ct/2, 
which contradicts f)5.3p for sufficiently large t. □ 

Appendix A. Construction of wave operators 

Here we briefly skecth a proof for the existence of the wave operators, or the 
solvability of the final state problem. For the construction of a nonlinear profile in 
the radial setting, we need only to consider a sequence of solutions in the form 

u n = U(t)f+f U{t - s)(nu,aD\u\ 2 )ds, (A.l) 

•J — t n 

with t n — > ±oo, which is normally transformed into such a form as 

u n =U(t)f- U(t + t n )B(U(-t n )f) + B(u n ) + Q- tn {un) + T_ tn (u n ), (A.2) 

where Q-t n and T_f n denote respectively Q and T with the Duhamel integration f* , 
and arbitrarily fixed (3, say (3 = 10. Below we consider only the case t n — > — oo, 
since the other case is similar. The following is the precise statement that we need 
for the nonlinear profile in this case. 
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Lemma A.l. Let f G H l x L 2 , R 3 t n — > — oo, and let {u n } be the sequence of 
solutions to the Zakharov system with the Cauchy data u n (—t n ) = U(—t n )f. Then 
there exist TeE and a unique u G S(T, oo) satisfying 

u = U(t)f + B{u) + Q^iu) + Too(m), lim \\u n - w|| S (r,oo) = 0, (A.3) 

as well as the Zakharov system on (T, oo). Moreover, if {u n } is bounded in L°°(M; H 1 x 
L 2 ), then u is global and the above convergence holds for any Tel. 

Proof. First, we can solve flA.31) on (T, oo) for T 3> 1, by the iteration argument 
similar to [7] in the space 

X := {u G C([T, oo); x L 2 ) | ||n|| s(Tj0o) < 1, ||iT||z(T,oo) < v}, (A.4) 

with ?7 := 2||[/(t)/||^( Too ) 1, using the estimates similar to (I4.54p as well as 

\\U(t)f\\ z(Ti00) < ||f/(t)/|| Lr>T (L6xH-) ^ as T ^ oo. (A.5) 

Similar estimates imply that u n are scattering as t — > oo for large n. Also similarly 
to (14 .54 p . we have for some 9 > 

\\U(t + t n )B(U(-t n )f)\\ s < \\B(U(-t n )f)\\ mxL * 

< \\u(-t n )n^ L 4u(-t n )fXe xH -^ o. ,6) 

Then by applying (14.541) to the difference equation, we obtain the convergence u n — > 
u in S(T,oo). Since t? n solves the Zakharov system, so does the limit u. If the 
former is uniformly bounded in H 1 x L 2 , so is the latter, and the convergence is also 
extended to arbitrary (T, oo) by the local wellposedness. □ 
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